Majorana zero modes bound to a vortex line in a topological superconductor 
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We explore Majorana zero modes bound to a vortex line in a three dimensional topological su- 
perconductor model, focusing our attention on the validity of the index theorem previously derived. 
We first solve the Bogoliubov-de Gennes equation at the zero energy to obtain the analytical index. 
We next calculate the topological index given by the order parameters. It turns out that they indeed 
coincide and that index theorem, which has been derived on the implicit assumption that a defect 
is point-like, is also valid for a line defect. 
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I. INTRODUCTION 

There has been renewed interest recently in topolog- 
ical phases of matter, owing to the rapid progress in 
both theoretical and experimental studies of topological 
insulators^ While the quantum Hall effect is a typical 
example of such a phase^i 3 - recent discovery is that even 
with unbroken time reversal symmetry, band insulators 
can be topologically nontriviali^— By extending the idea 
of the topological insulators to superconductors^ it has 
been shown that various topological phases are classified 
by ten universality classes, similarly to the random ma- 
trix theory ! 10 : 11 

Superconducting states are described by Bogoliubov- 
de Gennes Hamiltonians with particle-hole symmetry. 
This symmetry guarantees a pair of states with energy 
±E for every eigenstate. The zero energy states then 
have a special property which enables us to regard them 
as Majorana fermion modes. With many Majorana zero 
modes, the zero energy ground states are degenerate 
and these fermion modes are expected to obey the non- 
Abelian statistics. They are predicted in many other sys- 
tems such as a p-wave superconductor, etcj^^r— Generic 
conditions of the existence of the Majorana zero modes 
have been discussed i^£r— In particular, an index theorem 
has been derived ) 36 : 40 which reveals the relationship be- 
tween the index of the Majorana zero modes (analytical 
index) and a topological invariant given by the configu- 
ration of the order parameters (topological index) . 

Although previous studies have been focused on the 
non-Abelian statistics in lower spatial dimensions, Teo 
and Kane^ 3 - have recently proposed a three dimensional 
(3D) topological superconductor with a Majorana zero 
mode and discussed the possibility of the non-Abelian 
statistics in three dimensions. The index theorem has 
also been derived for the same system^ Since a point 
defect has been considered implicitly in these works, it 
may be interesting to consider a line defect inherent in 
three dimensions and to study whether Majorana modes 
appear or not in the core of a vortex line. It is also non- 
trivial whether the 3D topological index which is generi- 
cally the winding number of the order parameter on the 
two dimensional sphere gives a correct answer for a line 
defect which is expected to be classified by the winding 



number around a circle. 

In this paper, we study a 3D topological superconduc- 
tor model with a vortex line, focusing our attention on 
the validity of the index theorem derived previously when 
it is applied to the case of a line defect. In the next sec- 
tion, we discuss the model and its symmetry properties. 
In Sec. IIIIl we solve the Bogoliubov-de Gennes equa- 
tion analytically to show that in the core of the vortex 
line, there appear Majorana bound states exactly at the 
zero energy. We can then obtain the analytical index by 
counting the number of zero modes and their chirality. 
In Sec. IIV1 we calculate the topological index in three 
dimensions in the case of a line defect. We show that 
both the indices coincide indeed. Sec. [V] is devoted to 
summary and discussions. 



II. A MODEL AND ITS SYMMETRY 

We study in this paper a minimum model of a 3D 
topological superconductor proposed by Teo and Kane. 33 
This model includes a mass controlling the band inver- 
sion, which is assumed to change its sign along the z- 
axis. When the system is in the normal phase, an or- 
dinary insulating state appears in the positive mass do- 
main, whereas a topological insulating state is realized 
in the negative mass domain. At the boundary of these 
domains, massless Dirac fermions without doubling can 
appear. As shown by Teo and Kane, when such a system 
is in the superconducting phase, a topological defect such 
as a 't Hooft-Polyakov monopole background field yields 
a Majorana zero mode. In what follows, we consider not 
a point defect but a line defect, i.e., a vortex line, and 
show that the Majorana zero modes also appear in this 
model. 

To be concrete, we study the following Hamiltonian, 



rl — Hq -\-rLxi 



(2.1) 



where j = 1,2,3, and the 7 matrices obey {7%7 J } 
2S ij for i,j = 1,2, •••,6. We also define T 5 
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(— i) 3 7 1 7 2 ' ' ' 7 6 - To solve the Bogoliubov-de Gennes 
equation for zero modes, it may be convenient to choose 
the 7 matrices such that 

7 J = a ® a 3 ®1, 

7J+ 3 = cr 2 <g> 1 ® a\ 

r 5 = cr 3 ®l(g)l, (2.2) 

which is different from but unitary-equivalent to the rep- 
resentation used in Rcf.— . The term Hi in (|2.ip denotes 
the chemical potential term. 

To describe a vortex line, we introduce the cylindrical 
coordinates (r, 6, z). Let 0i and 02 be the set of the 
superconducting order parameter with a vortex which 
is a continuous function of r, 9, and z. We assume its 
asymptotic behavior as 

U «<, ^ _v / A ° ( cos q9 > sin ^ ( r ^ °°) O ^ 
\ A(r) (cos q6, sin q8) (\z\ 00) ' 1 Jj 

where the consistency at r — > 00 and |z| — » 00 requires 
that A(oo) = Ao- Without loss of generality, we assume 
that A > 0. Eq. ([273]) tells that at r -> 00, (0i,02) 
defines a map from S 1 to S . Therefore, for the function 
(0ij02) to be continuous, some zero points, (0i,02) = 
(0, 0), are indispensable somewhere in the (x, y) plane for 
every 2. In particular, in the limit \z\ — ¥ 00, we assume 
A(0) = 0. Therefore, these zeros form lines connecting 
r = at z — > ±00. As we shall show below, the details 
of these lines, for example, the shape of the lines and/or 
the number of the lines, are irrelevant to the index. A 
mass controlling the band inversion is, on the other hand, 
defined by 

03 = 03 (*)• (2.4) 

The mass is assumed asymptotically to be 

3 (±oo) = ±m . (2.5) 

As we shall see later, only the sign of mo is relevant to 
the index. In passing, we would like to mention that we 
have neglected, for simplicity, the vector potential of an 
external magnetic field which produces a vortex. 

Next, let us examine the symmetry of the model. De- 
fine the particle-hole conjugation operator 

C = i7 2 7 5 X, C 2 = 1. (2.6) 

Then, for the Hamiltonian (|2.1[) . we see 

CHC- 1 = -H, (2.7) 

which implies that the model belongs to the class D, a 
most generic class of superconductors.— ~— If the chemical 
potential Hi is neglected, the Hamiltonian Ho has higher 
symmetry. To see this, note 

T5H0T5 = —Ho- (2-8) 



This tells chiral symmetry of Ho, and the nonzero chem- 
ical potential term Hi breaks it. These twofold symme- 
tries of Ho lead to the following fictitious time reversal 
operator, 

T = r 5 c. 

From the transformation properties (|2.7p and (|2.8[) . as 
well as T§C = CT$ , it is readily seen that 

THqT- 1 = Ho, T 2 = 1. 

Therefore, the model without the chemical potential term 
Ho belongs to class BDl4i 

Although one may be interested in the Majorana zero 
modes mainly in the generic class of superconductors, i.e., 
in class D, the additional chiral symmetry (|2.8|) plays 
a crucial role in the discussion on the index theorem. 
Indeed, one can obtain zero modes even in the case of 
nonzero chemical potential, but without chiral symmetry, 
it is not possible to define the index any longer. There- 
fore, in the next section, we shall solve the zero modes 
separately for Ho and H. 

III. MAJORANA ZERO MODES 

In order to solve the Bogoliubov-de Gennes equation, 
we specify in this section the superconducting order pa- 
rameter by extrapolating the boundary value (|2.3j) to 
r = and/or z = simply as, 

(0 1 ,0 2 )=A(r)(cos^,sin^). (3.1) 

Namely, we assume that Eq. (|3.1|) is valid for any r, 
6, and z. We also assume that A(oo) = Ao, and that 
A(0) = 0, implying that the core of a vortex is located 
just on the z axis. In this case, the Hamiltonian is de- 
coupled into two parts, one depending on r and 9, and 
the other depending on z only. Such a model can be re- 
garded as a quasi 2D model interpolating the 2D model 
proposed by Jackiw and Rossi— and 3D model by Jackiw 
and Rebbi^ 3 - 

We also note that the order parameters (|3.1|) yields 
a continuous dynamical symmetry, which simplifies the 
calculations below. Namely, the order parameters Eqs. 
(|3.1j) and (|2.4p , and thus the Hamiltonian H are invariant 
under the rotation around the z axis. Indeed, we see 

[H,J 3 ]=0, (3.2) 

where the angular momentum around the z axis J3 is 
defined by 

J3 = -id e + - (1 <£> a 3 <£> 1 + ql ® 1 ® cr 3 ) 

= -i^-^(7V + 97 4 7 5 )- (3.3) 

Note also that [I^,^] = 0. It follows that the zero 
modes can be labeled by the chirality and the angular 
momentum. 
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In what follows, let us rewrite the Hamiltonian in a 
more convenient way to obtain the zero mode wave func- 
tions. The representation of the 7 matrices in Eq. (|2.2j) 
enables us to write the Hamiltonian (|2.1[) as 

( -fil ® cr 3 W- \ 
H \ \D+ -fil ® cr 3 J ' 

L> ± = a j ® Id, ± 1 ®cr%. 

Correspondingly, the particle-hole conjugation operator 
(|2.6|) can be denoted as 

C = icr 2 ®ia 2 K, (3.4) 

by means of which particle-hole symmetry (I2.7[) can be 
described such that 

CD ± C- 1 = -D ± . 

CltgJCT 3 ^ 1 = -1®(T 3 . 

The angular momentum J3 in Eq. (|3.3p is also decom- 
posed into 

^8 = 1® J3, 

where, J3 is the block diagonal element of the angular 
momentum operator defined by 

= -idg + S3, 
S3 = - (cr 3 ® 1 + ql ® cr 3 ) . 

53 is the sum of the real spin and the pseudo-spin due 
to the particle-hole space of the Nambu notation. The 
commutation relation (|3.2|) is translated into 

[D ± ,J 3 ] = \ M l®o a ,J 3 ]=0. (3.5) 

Now define 

U m (9) = e ^ m - s ^ e , (3.6) 



Di = u-\e)D ± u m {6) 

±A 

- p. 2m-(l+g) 
Wr 2r 

V 



Therefore, the transformed operator can be regarded 
as the operator with a definite angular momentum m. 



then, we see 

J 3 U m (6) = mU m (e). 

This matrix shall be used to construct the eigen func- 
tions of the zero modes. Therefore, the single- valuedness 
of U m (9) under the 2n rotation with respect to 9 requires 
that m should take in = (q + l)/2 mod 1. This implies 
that a unit vortex with q = 1 serves as a spin-1/2 ob- 
ject, and the total spin is just the sum of the spin and 
the pseudo-spin. The stability of the Majorana modes 
against the finite chemical potential depends on whether 
m is an integer or half-integer, as we shall see momentar- 
ily. It should also be noted that C and J 3 anticommute: 
CJ 3 = -J3C. It follows that 

CUmC- 1 = U- m . (3.7) 

Therefore, the zero modes appear as a pair of states with 
angular momentum ±m except for m = 0. This excep- 
tion plays a crucial role in the stability of the Majorana 
zero modes mentioned above. 

A. The case of fi — 

Since the model has chiral symmetry in the case of fi — 
0, the Bogoliubov-de Gennes equation at the zero energy 
is decoupled into two sectors with definite chiralities. Set 
the wave function as * m = '^!^ n ) T , where ± denote 
the chirality of the zero mode wave functions. Then, the 
equation for the zero mode can be written as 



Because of Eq. (|3.5[) , the zero modes can be simultaneous 
eigenstates of the chirality and the angular momentum 
J 3 . Therefore, we have introduced the quantum number 
m which labels the angular momentum, Jz^!^ — wity^. 
It is easy to see that the unitary transformation U m {9) 
defined in Eq. (|3. 61) removes the ^-dependence from : 



° \ 

2m+(l+q) 
±A 

I 

Correspondingly, the eigenstate of J3 is defined as 
*±(r,9,z) = U m (8)^Jr,z). 



r 



±A 

^3 T 03 


2m-(l-g) 
2r 



p, 2m+(l-g) 



~(d 3 T0s) 
±A 
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Thus, it turns out that the ^-dependence of the zero mode 
wave functions is determined by U m . 

Next, note that the z-dependence is only from the di- 
agonal elements of the Hamiltonian. To solve it, we set 



Vl\ 



Ill * Vim' V2m' ^2m) ■ 



Then, we have 

where we have used the notation Vm( r ) = 
(vL(r),vL(r)) T and £±(r) - (tL(r),& m (r)) T ■ 
The asymptotic mass to in Eq. (|2.5I) tells that the 
above two types of the wave functions, r/ and £, are 
normalizable only for mo > and mo < 0, respectively. 
Thus, the normalizability of the wave functions in the 
limit \z\ oo requires that the zero mode equation 
should reduce to 



of the two components rjf m and r)f m for each allowed 
m. This parameter plays a role in the normalizability 
in the limit, r — > oo. In this limit, 1/r terms in Eq. 
(|3.8p can be neglected, and it is easy to show that the 
solution of the equation is a linear combination of e ±A ° r . 
Therefore, one can fine-tune the remaining parameter to 
vanish the exponentially growing component e +A ° r . This 
is the normalizable wave function at r — > oo. 

Final results, including the case of too < 0, are sum- 
marized in Tabic [T] 



B. The case of fx 

Since chiral symmetry is broken by the nonzero chemi- 
cal potential, the wave functions ijj^ (r, z) couple together 
when /i 7^ 0. Even in this case, however, rj and £ do not 
coexist in the same wave function, owing to the normal- 
izability in the limit \z\ — > oo. The zero mode equation 
is thus given by 



£(»•)= 
V 

where 



± (r) = 0, 



± ; 



dt ~£±(r) = 



(to > 0), 
(to < 0), 



(3.8) 



d ± = 



Or 




Finally, we examine the normalizability in the limit 
r — » and oo. Let us consider the case of m > 0, for 
example. In the vicinity of r = 0, we can set A(r) ~ 0, 
which implies that rf^ m and r^ m are decoupled in Eq. 
(|3.8p . Then, we see that the wave function behaves such 
that r)f m ~ r ™-0-±i)/* and r)f m ~ r -™-(i±<z)/2. F or the 
wave function to be normalizable, we should impose the 
condition that ^l^^l 2 oc r a with a > — 1 for both j = 
1, 2. We then see that 77^ are normalizable, respectively, 
only for \m\ < (1 =p <?)/2. Therefore, in the case of q > 0, 
we should impose rj m — 0, since 77" is normalizable but 
77^ is unnormalizable. Likewise, in the case of q < 0, we 
should impose 77" = 0. 



m 


q 


Nonzero wave fn a 


Index 6 


+ 


+ 




-9 




+ 


£+ 

S> 771 


+q 


+ 




Vm 


\q\ 








~\q\ 



"The angular momentum m takes m 



kl-i 



2ii mod 1 in 



< 



. Therefore, in each case, there are \q\ normalizable zero 



m < 
modes. 

'The index of the Hamiltonian is defined in Eq. 1 14.111 

TABLE I: Nonzero components of the normalizable wave 
functions 

It should be noted that the wave function discussed 
above has still an ambiguity with respect to the ratio 



0. 



Vm = °> 



d + 

d 



ificr 





( Vm 


—i/ia 3 ) 


V Vm 






i/jfj 3 J 





= 0, (to > 0), 

= 0, (too < 0). 

(3.9) 



These equations do not allow normalizable solutions in 
general. However, one exception is the case m = 0, 
which occurs only when q — odd. To see this, note that 
Eqs. (|3.4p and (|3.7p enable us to impose the relationship 
±C'0m( r j z ) = Tp-mi 7 "' z )- Hence, m = is quite special 
in the sense that independent variables (the components 
of rj^ or are reduced: rjf = 77^* and £^ 
Therefore, we have for too > 

1 + <? N 



"sio • 



r 
8 r 

and for too < 



2r 
2r 



V10 + A V 



10 



Vw - A77. 



10 



W10 



+* 
W10 



1 + g 

2r 

1-g 
2r 



er +Aer *- 



ei + o-ACib* 



^£10* 



0. 

= 0, 

0, 

0. 



The normalizable solution is for toq > 



Vw( r ) 
Vw(r) 



and for toq < 



oc e 



- f drA(r) 



Jl 



-iJ(g_i)/2(M?") / 



(9+l)/2 



oc e 



-/ r drA(r) / J(q+l)/2 



(?-l)/2 



( M r) i ' 
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where J„(x) is the Bessel function of the first kind. This 
result implies that if one adds chiral symmetry breaking 
perturbations, only one zero mode survives when q is odd, 
whereas no zero modes survive when q is even. In this 
sense, the Majorana zero mode in class D is classified by 
Z 2 16 - 17 - 36 - 44 

In passing, we mention that the zero mode solution ob- 
tained in this subsection is nothing to do with the index, 
since it cannot be an eigen state of the chirality T5. 

IV. INDEX THEOREM 

The index is associated with the zero modes with def- 
inite chiralities, reflecting an analytical property of the 
Hamiltonian as a differential operator, and the index the- 
orem relates the index with a topological invariant given 
by the order parameter. In this section, we first review 
the index theorem obtained previously, 3 - and next show 
that the analytical index and the topological index of the 
present model indeed coincide exactly. 

It should be noted that chiral symmetry plays a central 
role in the index theorem. Therefore, we restrict our 
discussions in this section to the Hamiltonian 

A. Summary of index theorem 

The index of Ho is defined as 

indWo=JV+-JV_, (4.1) 

where N± is respectively the number of zero modes with 
chirality ±. The chirality is here defined as the eigenvalue 
of IV From the zero modes obtained in Sec. IIII1 we can 
compute the index, which is summarized in Table HI It 
has been shown that the index of Hq can be expressed 

indfto = ~ J dSjJ j {x,0,oo), (4.2) 

where dSj is the infinitesimal surface element at the 
boundary of R d (typically S^" 1 at r — > oo), and 
Ji(x, Mq, Moo) is the axial- vector current defined by 

J 3 "(aj,M ,M oo ) 

= lim tacyRV ( | 8{x — y). 

y^x 15 ' \Mo-iH Moo-mJ v y) 

Since the current itself is not well-defined, a Pauli-Villars 
regulator with a mass M^ has been taken into account. 
In three dimensions, we have ) 36 i 46 in the limits Mo — > 

and Mqo — > °o 

J l (x,0,oo) = J—e^e^adjfodkfa, 



where i,..,a,.. = 1,2,3 and if) 2 — Y^j4>]- In Refj2£, we 
assumed that if) 2 is constant at t he boundary S 2 of R 3 
(in the limit r = \J x 2 + y 2 + z 2 — > oo), which implies 
that the singularity is due to a point defect and the sys- 
tem is spherically symmetric at infinity. In the present 
case, however, the vortex is a line defect in the super- 
conducting order parameter {<f>\, <f>2 )j and the system has 
cylindrical symmetry rather than the spherical symme- 
try. Therefore, 4> 2 cannot be uniform at the boundary 
r = \J x 2 + y 2 — y oo or | z | — y oo . The index theorem 
(I4.2[) is nevertheless valid if the integral is performed over 
an appropriate surface, a cylinder for the present model, 
reflecting the symmetry of the model. 



B. Calculation of topological index 

It is noted that the topological index, the right hand 
side of Eq. (|4.2[) , is determined only by the bound- 
ary value of the order parameter. Therefore, in the fol- 
lowing discussions, the superconducting order parameter 
(<^i)<^2) is n °t necessarily given by Eq. (|3.1[) : Only the 
assumption of Eq. (12.31) is enough to compute the in- 
dex, and therefore, even if the a vortex line is curved 
or deformed, the topological index is the same, and so 
should be the analytical index, provided that the order 
parameter obeys Eq. (|2.3|) at the boundary. This is one 
of advantages to utilize the topological index rather that 
the analytic index. 

In order to calculate the integral over a cylindri- 
cal surface in Eq. (I4.2[) . it should be noted that in 
terms of the form, the surface element dSj is given 
by the following two-form, dSj = (l/2\)ejkidx k dx l , 
where dx k dx l — —dx l dx k , and that the order param- 
eter depends on the cylindrical coordinates such that 
4> = (cj)i(r,8) } (f>2(r,9) 7 (f}3(z)) at the boundary. On the 
(9, z) surface in the limit r — >■ oo, we have J J dSj = 
(J 1 da; 2 — J 2 dx x )dx 3 . Here, the current at the bound- 
ary r — > oo is given by 

3 = MAg + *§(z))3/» e e 

for i,j,a,b — 1,2, where <fr' 3 (z) = d z 4>z(z). On the other 
hand, on the (r, 9) surface in the limit z — > ±oo, we have 
J J dSj = J 3 dx 1 dx 2 , where 

3 -27(A 2 (r)+m 2 ) 3 / 2e M" d *'h» 
for a, b = 1, 2. Substituting these into Eq. (|4.2|) , we have 



6 



ind'Ho 



1 

-in 



dz 



(A§ + $(*)) 3 / a 



2 

47T 



d 2 x 



(A 2 (r) +to 2 ) 3 / 2 



e afc c>i</> a d 2 b , 



(4.3) 



where i,a, b = 1,2 and d 2 x = dx 1 dx 2 . The factor 2 in 
the latter term is due to contribution from the two limits, 
z — > ±oo. The integration over z in the former term can 
be carried out as follows: 



dz 



<j>'{z) 



(A 2 +02( 2 ))3/2 



dx(x 2 + A*)" 3 / 2 



2ron 



Ag^/mg + A 2 

Therefore, the first term becomes, after the integration 
over z. 



m 



1 



V ^+A|2 7 r 



dxV^S^f,, 



(4.4) 



where <j> = </>/Ao at r — oo. The second term in (|4.3p 
can be calculated as 



1 

2^ 



d0 



dr- 



m 



„ (A 2 (r) + to 2 ,) 3 / 2 



1 /^fl «*U « X r a to AA' 
2tt / J (A 2 (r) +mg) 3 / 2 

where we have used the relation <fi a (r,9) = 4> a (8)A(r), 
which is valid from Eq. (12.31) . The integral over r thus 
reduces to 



dr- 



m AA' 



dir- 



ra 



o (A 2 (r)+TO 2 ) 3 / 2 7 2(x + m 2 )3/ 2 

m m 



V™o \/Aq + r, 
Therefore, the second term in Eq. (|4.3j) becomes 



sgn(m ) 



rn 



2vr 



d0 e a Va9^fc. (4.5) 



From Eqs. fO} and ([33]) . it turns out that Eq. P~5|) 
finally becomes 



ind'Ho = -sgn(mo)— rf> dde ab § a de4> h . 

Z7T 

= -sgn(m )g. 



(4.6) 



The right hand side of Eq. (|4.6|) , the topological index, 
is indeed the same as the analytical index, the left hand 
side, summarized in Table HI We also note that the first 
line of Eq. (|4.6|) tells that the present index is a product 
of the sign of m and the winding number around a vortex 
which is a topological invariant of a line defect. 



V. SUMMARY AND DISCUSSIONS 

We have studied a 3D model of a topological super- 
conductor with a vortex line. In addition to particle- 
hole symmetry which is the most basic symmetry of the 
Bogoliubov-de Gennes Hamiltonian, the model has chi- 
ral symmetry when the chemical potential is zero. This 
enables us to explore the Majorana zero modes in terms 
of the index theorem. 

In the former part of the paper, we have solved the 
Bogoliubov-de Gennes equation at the zero energy for 
the Hamiltonian with a vorticity-g vortex. When the 
chemical potential is zero, we have obtained \q\ Majorana 
zero modes. 

To solve the Bogoliubov-de Gennes equation, we have 
made the best use of the conserved angular momentum 
around the z axis. Precisely speaking, it may not be 
suitable to call the zero modes thus obtained Majorana 
modes, since it is not "real" . To see this, note that Eq. 
(I3.7[) allows us to choose ±C\E'^ = *^ m . Namely, "com- 
plex conjugate" of VP^ is not itself but is which 
is orthogonal to In this case, linear combinations 

(*£ + *iJ/V2 and - *tJ/V2, and similarly 

for ^± m can be regarded as Majorana modes. Alterna- 
tively, one can add rotational symmetry breaking per- 
turbations. Then, to is not a good quantum number any 
longer, and one can obtain generically Majorana wave 
functions. Even in such a case, the index theorem guar- 
antees that the Majorana zero modes are topologically 
protected and such small perturbations cannot change 
the index, which implies at least \q\ zero modes even for 
a perturbed model. 

Using \q\ zero modes obtained so far and taking their 
chirality into account, we have calculated the analytical 
index. On the other hand, in the case of nonzero chemi- 
cal potential, we have obtained just one zero mode when 
q is odd, whereas no zero modes when q is even. There- 
fore, this result suggests that Majorana zero modes are 
classified by Z for systems of class BDI, whereas by Z 2 
for systems of class D i 16 i 17 ' 36 i 44 

This classification scheme can also be supported by the 
following arguments 4£ Let us consider the case too > 
and q > 0, for example. When [i = 0, we have q 
zero modes whose eigen functions are denoted as ^ m = 
(0, \I/~) T , where \m\ < (q — l)/2. For these unperturbed 
states, we can compute the first order perturbative cor- 
rections of the chemical potential term, or more generi- 
cally, of some hermitian operator O which is odd under 
particle- hole transformation COC^ 1 = — O but even un- 
der chiral transformation — O. Note that we can 
choose the phase of the states such that C^ m = \&_ TO . 
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Therefore, we have (* m |0|* m ) = -(*_ m |C|*_ m ). 
These matrix elements do not vanish in general for oper- 
ators breaking chiral symmetry and we have ±e energy 
corrections to these states4£ However, there is one ex- 
ception m = which occurs when q is odd. The above 
equation readily leads to (^olCl^o) = 0. Therefore, we 
can expect that state is protected from chiral symme- 
try breaking perturbations. 

In the latter part of the paper, we have calculated the 
topological index. In the previous paper^ we have de- 
rived generic index theorem on the implicit assumption 
that the defect of the model is point-like. However, in 
3D superconductors, a vortex is a line defect. It has been 
shown that the index theorem obtained in the previous 
paper is valid even in the present case of a line defect. 

It turns out that the index theorem is a quite pow- 



erful tool to explore the Majorana zero modes bound 
to topological defects in superconductors. However, it 
needs chiral symmetry by definition. Therefore, it may 
be interesting to extend the index theorem such that it 
can be applicable to more generic universality classes of 
superconductors . 
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